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Abstract 



The behaviour of the flat anisotropic model of the Universe with a scalar field is 
explored within the framework of quantum cosmology. The principal moment of the 
account of an anisotropy is the presence either negative potential barrier or positive 
repelling wall. In the first case occur the above barrier refiection of the wave function 
of the Universe, in the second one there is bounce off a potential wall. The further 
evolution of the Universe represents an exponential inflating with fast losses of an 
anisotropy and approach to the standard cosmological scenario. 



1. Introduction 



One of the basic problems of the standard cosmological scenario is the presence of an 
initial singularity. For its elimination are used the various approaches, one of which is 
a quantum cosmology originated by DeWitt [1] more than thirty years ago. In last one 
the whole Universe is described with the help of the wave function. Until recently main 
part of papers was devoted to the consideration of the process of quantum creation 
of the closed Universe. In this direction by number of authors [2] was developed the 
scenario of the spontaneous creation of the closed Universe from "nothing", where 
"nothing" means the state with absence not only matter, but also the spacetime in 
classical understanding. Thus it is necessary to note that the closed Universe has a 
zero total energy and charge therefore there is no violation of any conservation laws. 
At describing of process of creation of the Universe the so-called instanton method [3] 
is widely used. In the last one the closure of the Universe is a necessary requirement 
of the quantum creation, otherwise the action is equal to infinity and, hence, relevant 
probability of the creation is near to zero. However, recently observant data have 
appeared according to which our Universe probably is open. A number of models of 
the quantum nucleating of the open Universe were proposed in this connection. One 
of its is the Hawking- Turok instanton [4], for which the opportunity of an analytic 
continuation of the instanton solution not only in the direction of creation of the 
closed Universe, but also in the direction of the open Universe is shown. This result, 
being interesting on itself, is not unexpected. For a long time it is known (see e.g. [5]) 
that dividing by various ways the four-dimensional de Sitter Universe on space and 
time we can obtain the various de Sitter Universes: closed, fiat or open. The quantum 
creation of the open Universe can seem not actual in view of the noted above infinite 
value of action. But here it is possible to remember the theory of nuclcation of the 
open Universe from bubble [6,7] filled by false vacuum, when the Universe inside of 
this bubble looks finite from the point of view of an external observer and open, that 
is infinite, from the point of view of an inside observer. Therefore, for an external 
observer the process of the instanton nucleation of the open Universe is quite possible 
and the growth of the bubble size is carried out by transferring of energy from the 
surrounding de Sitter space. 

2. Basic classical equations 

We shall consider an anisotropic Bianchi type I model. For it the synchronous form 
of metric can be expressed as (velocity of light is equal 1): 

ds'^ ^ - al{t) dx^ - al{t) dy'^ -al{t)dz^, ( 1 ) 
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where by ai,a2 and 03 the scale factors on directions x, and z accordingly are 
designated. This model is an anisotropic generalization of the Friedmann model with 
Euclidean spatial geometry. Three scale factors ai,a2 and 03 are determined via 
Einstein's equations. For convenience of realization of the analytical calculations we 
can write them as follows [8] : 

a\ = r{t)qi, aa = r(t)g2, cts = rit)q3 , (2) 

where gi, g2, are dimensionless variable subordinated to the following requirements: 



n 9" = 1' J2 i<ia/qa 



0, 



(3) 



a=l 



(dot means derivative with respect to time t); whence follows that tl = r^- 

For the line element (|1|) with the account of the components of the Ricci tensor 
write as: 
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(4) 



Let's use the last one for obtaining (1-1) and (2-2) components of the Einstein tensor: 
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Subtracting from G\ the component G2 one obtains: 
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Entering in the last equation the notification Qa/B = {q2/q2 — qi/qi) let's have: 



3^ + ^ = 0, 

r Qal3 
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that after integration gives: 
where Cap are integration constants. From here we get: 



Qa/3 — Cafs/ r^, 



Qa _ Co 



(5) 



and according to requirements (||), I] = 0. Thus integrating the last equation 

0=1 

one finds: 



Qa = A^exp J ^1 



(6) 



where Aa are integration constants and fi 

a=l 

from P) is got: 



1. Now, using the relation 
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a=l 



Co determines an anisotropy of the given model. 



a=l 



3. Quantum evolution of the Universe 

The possibility of the nucleation of an open Universe circumscribed in Introduction 
gives the basis for the consideration of the quantum creation of a fiat Universe. As is 
known the basic equation of the quantum cosmology is the Wheeler-DeWitt (WDW) 
equation. For its making we shall consider the theory of a scalar field ip with La- 
grangian 

L = -R/2 + {d^iff/2~Viip) (8) 
or using the expression for scalar curvature R from (H) one obtains 



„3 3 
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(9) 
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(accurate to complete derivative with respect to t). The relevant conjugate momen- 
tums are equal 



dL o dL 

0(p or 



and the Hamiltonian of the system is 



3 

H = p^ip+ Pri- + J2 PqJc. - (11) 
a=l 

Let's note, that for deriving the exact equations it was necessary to use the expression 
(I) for scalar curvature instead of (0). Using the last one is impossible whereas in 
this one the integration for elimination of (ja/Qa is already yielded. It is intolerable as 
actually there is the deletion of variables Qa and thus truncation of the Hamiltonian. 
Let's note also, that if to use the last of relations (^) and expression (ja/la = Ca/r^ 
from (^ exchanging simultaneously Pg^ (here Pg^ = —id/dqa), it is easy to 

obtain that QaPqJ^ = Co^- The last means that in our case is the eigenfunction 
of the operator pg^ . It allows with account (^) to write the Hamiltonian (|11]) as: 

^ = - ^ + r'V{^) + (E Cl) /2r\ (12) 

Quantizing (|1^) by replacement of momentums p^ and pr on —id /dip and —id /Or 
accordingly and also using the rescaling ip -\/6$ we obtain the Klein-Gordon equa- 
tion 

' 1 d ( d\ 1 92 1 

U^f = -^(jZCl^ It' - I2r'vm. (13) 

It is the required WDW equation in minisuperspace of the variables r and $. In this 
equation the parameter p represents the ambiguity in the ordering of noncommuting 
operators r and pr- Let's emphasize, that the wave function of the Universe ^ does 
not depend on time. This circumstance is valid for the closed Universe by virtue of 
equality to zero of its total energy remains valid and for the flat Universe on the basis 
of given in Introduction reasonings. 
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The transformation of the wave function 



allows to eliminate first derivative in (T^) 

1 



j,2 Q(^2 



u, 



ilj{r, $) = 0, 



P 



1 QECl 

- ^4 12rV($). 



(14) 



Let's note, that in expression for an effective potential energy there is an addend of 
"centrifugal energy" — p (1 — p/2) /2r^. 

Is explored Eq. (|1^) with a various form of the potential l^($) and by choice of 
parameter p. 

A. Above barrier reflection of the wave function of the Universe 

1. de Sitter minisuperspace. Let's consider the simplest case of minisuperspace model, 
when the factor ordering p = and potential V{^) represents constant vacuum energy 
creating an effective cosmo logical constant. Then WDW Eq. ([T^) will take the 
form of the one- dimensional Schrodinger equation 



ef 



^(r) = 0, 



(15) 



where = 12ey is the Hubble parameter. Entering p = H^'^r and 7 = 6 C*^ we 
can rewrite (|r^) as 



a=l 



dp^ 



ef 



^{p)=0, Uef = -l/p' 



(16) 



Eq. (|T6D describes the motion of a "particle" with zero-point energy in the field 
of the effective potential Uef- The interesting feature of the given potential is that 
near to the origin of coordinates it is approach infinity under the law Uef ~ —7 / p^ 
(that is we can neglect the second term in potential). As is known [9] , this case 
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is intermediate between when there are usual stationary states and the cases when 
happens the "collapse" of a particle in the origin of coordinates. Therefore, it is 
necessary to carry out the additional analysis here. 

For this purpose we shall search near p = for solution of 



as ~ that at substitution in gives 



*(P) = (17) 



— s + 7 



with roots 



51 = 1/2 + ^1/4-7, 52 = 1/2-^1/4-7- 



Further, we select around the origin of coordinates small area of radius po and is 
exchanged in it the function — 7/p^ by the constant — 7/Po- Having defined wave 
functions in such "cut off" field then we shall look what happens at the passage to 
the limit po — 0. 

Let's assume at first that 7 < 1/4. Then 

si > S2 > and at p > po the general solution of Eq. (|l^ looks like (at 
small p) 

m = Ap'' + Bp'^ (18) 
( A, B are the constants). At p < po the solution of 

d''^ 7 T n 
h — \E' = 0, 

dp^ p^ ' 

finiteness in the origin of coordinates looks like 

^ = Csin(A;p), k = ^ / p^. 

At p = Po function \& and its derivative should be continuous functions. Therefore, we 
can write one of requirements as the requirement of the continuity of the logarithmic 
derivative with respect to \E' that gives in 

s,pr'' + {B/A)s, 

^^^^= pr-^ + iB/A) ' 

or solving rather (B/A) we obtain 
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B / A = const pI^-'\ (19) 

Passing on to the limit po ^ we find that B /A ^ 0. Thus from two solutions (|18D 
remains 

^ = Ap'\ (20) 
Let now 7 > 1/4. Then Si and S2 are complex: 

si = -1/2 + 2^^7-1/4, S2 = sl. 

By analogy to the previous reasonings we again come to equality ([T9|) which at sub- 
stitution of values si and S2 gives 

B/A = const fl/^. (21) 

At po this expression does not approach to any limit so the direct passage to 
the limit po ^ is impossible. The general form of the real solution of (|l^) can be 
represented as following: 



= const i/pcos [^1 — 1/4 In (p/po) + const j . (22) 

This function has zeros which number unlimited grows with decreasing of po. Then 
at any finite value of the total energy E the "normal state" of a "particle" in given 
field corresponds to the energy E = —00. As "particle" is in infinitesimal area around 
of the origin of coordinates there is the "collapse" of the "particle" on centre. 

Further, we find the vector of the probability density flux near to zero. In our 
one-dimensional case we have: 

2\ dp dp )' 

It is easy to obtain from here that: 1) in the case 7 < 1/4 at use from ( pO] ) relevant 
probability density flux j = (as well as for any real wave function). 2) At 7 > 1/4 
with the account from ( ^2]) we have j = — 1/4. The upper sign corresponds 
to the ingoing and the lower one to the outgoing wave. The obtained result means 
that there is the constant probability density flux near to the origin of coordinates. 

Thus we have two types of the behaviour of the wave function of a "particle" at 
various values of the parameter of an anisotropy 7: 1) at 7 < 1/4 the wave function 
near to the origin of coordinates tent to zero; 2) at 7 > 1/4 happens the "collapse" 
of a "particle" on centre. 
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Figure 1: The above harrier reflection of wave function: 1 are ingoing and reflected 
waves; 2 is the passed wave. 



Let's note that the probleme of falhng of a particle on centre has been aheady 
studied in isotropic model with the matter equation of state p = e in [10]. Let's con- 
sider further the motion of a "particle" in the semiclassical approximation. From the 
form of potential (|T6D and equality to zero of the total energy of the Universe follows 
that the examination of Eq. (|16D is reduced to the one-dimensional problem about 
above barrier reflection (Fig. 1), i.e. to reflection of a "particle" with energy exceed- 
ing height of the barrier. In our case "particle" goes down to some "turning point" p2 
in which it changes the direction of motion on an inverse. The given point represents 
a complex solution of the equation Uef = 0, namely p2 = 7^/^ exp (?7r(2n + l)/6). 
Then the required reflectance R is found as [9]: 



R = exp[ -4 Im / [-U, 



P2 



Pi 



'ef\ 



,1/2 



(23) 



(here pi is any point on the real axis). In the considered case using ( |T6| ) and (^) we 
have: 

i? = expf-^V7vrV (24) 



The last expression is obtained with use of the semiclassical approximation. It is 



useful also to find the exact solution of Eq. (|l^) which looks like: 

*(P) = Civ/pJ (^^1/4- 7, \p') + C^^pY (^^1/4- 7, \p') , (25) 

where Ci and C2 are integration constants, J(z/, z), y(z/, z) - the Bessel functions of the 
first and second kind accordingly, z/, z are index and argument of these functions. For 
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the finding of tlie reflectance it is necessary to consider the behaviour of this solution 
on major distances from the origin of coordinates. In this case ^'(p) describes ingoing 
and reflected waves and the reflectance will be equal to the ratio of amplitudes of 
these waves. The asymptotic form of the function \E'(p) at p — *• oo is 

(26) 
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for 7 < 1/4 : R 



It is clear from here that: 



for 7 > 1/4 : R = exp (^--7r^7 - l/4j . (27) 

The second expression from (^) at 7 ^ 1/4 coincides with reflectance from the 
semiclassical approximation (l2^). 

Being based on results obtained at the analysis of behaviour of the wave function 
of the Universe near to the origin of coordinates is concluded that at 7 < 1/4 does 
not happen the accumulation of the wave function at p ^ and consequently takes 
place the complete reflection of the wave function from the barrier {R = 1). In the 
case 7 > 1/4 which corresponds to collapse of a "particle" on centre the reflectance 
R becomes less than 1. It happens because there is nonzero probability density flux 
in the infinitesimal area around of the origin of coordinates. Let's note that at the 
approach to zero the problem, generally speaking, ceases to be stationary. It gives 
that the wave function can accumulate in this area. 



2. Variable scalar field. Let's consider further the case of Eq. (|14|) when the factor 
ordering p is equal to zero as before, but potential of the scalar field is variable: 
V{(p) = Introducing the rescaling ip -\/6$ and m ^ p/6 we shall write 

the WDW equation as: 



^ = 0, 



f/e/ = -J-pV$l (28) 

The finding of the analytical solution of the obtained equation represents a complex 
problem. Therefore, for its examination we shall search the WKB solution as '^c = 
e~^^ . The relevant equation for action S{r, $) is 

fdsV 1 fdsV , , 
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For finding of tlie solution of this nonlinear differential equation it is possible to reduce 
it to system of the ordinary differential equations called the characteristic system of 
the given partial equation. Using this system it is possible to construct an integrated 



surface of Eq. (^91) , consisting from the characteristics. The required system of the 



characteristics written with respect to arbitrary parameter t has a form [11] 

r\t) = -2p, $'(t) = lg, 



S'it) = -2 4 + A'<^' 



Pit) = ^g' + 2^ - Vr3<I)2, (30) 
q{t) = -2/iV^$, 



where ' denote a derivative with respect to t and introducing denotations p = dS/dr, 
q = dS/d^. The obtained system of equations describes an one- dimensional motion 
of a "particle" along the characteristic. In this case monotonically varying parameter 
t can play a role of time [12] in due course there is an evolution of the Universe. 

The obtained system of characteristics can be use for finding of dependence of 
coefficient of above barrier reflection R of wave function of the Universe from value 
of the field 

But before we shall note one useful possibility of simplification of making of a 
calculation. Expand an effective potential Uef from (|16D near to a maximum in a 
Taylor series: 

f^e/ = f^max + - (P - Pmax)^ (31) 

here f/max and pmax are a maximal value of potential and relevant value of p, and a 
is a value of a second derivative with respect to Uef from (|1^) in the point pmax- The 
values of specified quantities are easily finding from 



a = -24 



Then using (^3]) and (|3T|) we have 



^ = exp(-y^v^7r|, (32) 
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Figure 2: Dependence of coefficient of above barrier reflection R from value of the 
scalar field $. 



that approximately coincides with (p^ . Thus in an one-dimensional case there is an 
possibility to use approximate expression (pT]) instead of exact one. 

Further, we shall search for the solution of system (|30|) numerically. Using ( p3|) 
and (pTD it is possible to obtain the following form of dependence -R($) (Fig. 2). It is 
obvious that the reflectivity promptly decreases with increasing of value of the field 
$. After reflection of the wave function the following necessary stage of evolution of 
the Universe is the inflationary period providing the "stretch" of linear size of the 
Universe with Planck up to macroscopic. For ensuring of the sufficiently long period 
of inflation it is necessary the realization of two requirements on scalar field (see e.g. 
[13] ): 1) it should be a Planck order; 2) its values should vary slowly with time. 
The realization first of these requirements corresponds rather small, but nevertheless 
distinct from zero refiectivity R (see Fig. 2). Solving system (^) will be easily 
convinced that is fulfilled second of the posed requirements also. 

B. The wave function bouncing off a potential wall 

Now we shall consider Eq. (|n|) in a case when the parameter p is different from 
zero and there is a potential of the scalar field V{(p) = At such setting 

of a problem we have two essentially various variant of the effective potential f/e/- 
At realization of a requirement < p < 2 the form of Uef from ([M]) as a matter 
of fact by nothing differs from the case when p = and the relevant consideration 
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4 2 

Figure 3: The effective potential Uef ■ I is a classically forbidden area, II is a classically 
allowed area. 

will be made by analogy to the problem studied in the previous section. In a case 
omissions of the mentioned above requirement the qualitatively new statement of a 
problem is possible. For this purpose the realization of one requirement is necessary: 
|— + p/2\ > 7, which ensures occurrence in the effective potential of a positive 
factor before greater than 7. Thus there is a possibility of occurrence of a 
repelling potential wall. Then the effective potential in (|l^) will be: 

£/.; = tP!Z^_^_,V*^_, ,33) 

Thus the influence of a massless scalar field (e. g. photon gas) with energy density 
e also is taken into account, the sense of which introduction will be explained below. 
The form of potential is shown on Fig. 3. Then Eq. (|l^) describes a motion of a 
"particle" with a zero-point energy in potential Uef{r, $). 

As is obvious from Fig. 3 the effective potential parts space on two areas: classically 
forbidden (interior) and classically allowed (exterior). At transition of a "particle" 
in the I-st area its wave function will be exponentially decay (as Uef{r, $) approach 
infinity at r ^ 00) and the relevant probability of realization of the given state tends 
to zero. 

On the other hand, the wave function \1/ describes the wave incident on the barrier 
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Uef{r, $) on the part of major r and the wave reflex from the barrier. The physical 
sense of the given statement consists that the incident wave describes the contracting 
Universe and reflex - expanding one. 

We shall consider an evolution of the scalar field $ at stages of contraction and 




expansion. For this purpose we shall search the WKB solution of the equation (|T^ 
with potential (|33|) in classically allowed range as \E'c = e"*"^. The relevant equation 
for action S{r, $) will be: 

= 0. (34) 

By analogy to the case of Eq. (|29|) we shall make system of the characteristics with 
respect to arbitrary parameter t: 

r'{t) = -2p, <^'{t) = ^q, 



p'(t) = P/2 + P/2| ^2^-4/iW, (35) 



As well as in case of Eq. (p9D system of equations (|35D is describes an one-dimensional 
motion of a "particle" along the characteristic. Thus the behaviour of a "particle" is 
similar on bounce off a potential barrier $(r) made at cross of the effective potential 
with the plane $ — r (that is when Uef = 0). Therefore, a semiclassical wave function 
\E'c describes an ensemble of classical universes evolving along the characteristics S. 
Then the ensemble of these characteristics can be considered as the trajectories of a 
motion with the various initial conditions. 

We shall note here that introduced earlier the massless scalar field with the energy 
density e is necessary that the Universe having a zero total energy always remained 
in classically allowed area (Fig. 4). 

As it was specified above the evolution of the Universe is described by two stages: 
at first regime of contracting and then, after bounce off, regime of expansion. At 
the stage of contracting the field makes oscillations with increasing amplitude. After 
bounce off two variants are possible: 1) the field $ amount to rather major value and 
after reflection varies slowly that corresponds to an inflationary period and further 
transfers on scalaron stage (Fig. 5); 2) $ hasn't amount to major values and having 
reflected at once makes fast oscillations losing the energy. Thus, naturally, there is 
no inflationary stage (or it is too short). 
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Figure 4. Conventional 
form of the potential barrier 
with 6. I is a classically 
forbidden area; II is a clas- 
sically allowed area. The 
dashed lines - form of the 
barrier without e. 



Figure 5. The diagram of dependence of the scalar 
field $ from r. 1 is an evolution of the field on 
the contracting stage, 2 - during the inflation. The 
dashed line conventionally shows the form of the po- 
tential barrier. 



It is necessary apart to note that in considered model the presence of a repelhng 
potential wall does not allow the Universe to collapse. 



4. Conclusions 



In the submitted paper we have considered an anisotropic cosmological Bianchi type 
I model with the scalar field. Distinctive feature at the solution of the given problem 
was representation of three scale factors in such manner that the final Einstein's 
equations have turned out dependent only from one function r(t). 

In the quantum approach the basic equation of the quantum cosmology - WDW 
equation (|T^) was obtained. The examination of the latter was reduced to considering 
the following problems: 

1. The simplest model with constant scalar field playing a role of an effective 
cosmological constant was studied. It has allowed to reduce the WDW equation 
to the one-dimensional Schrodinger equation. The presence of the parameter of an 
anisotropy 7 in the effective potential gives in interesting feature: there is some critical 
value of this parameter at which there is partitioning the problem into two variants. 
In the first case at 7 > 1/4 the collapse of the wave function on centre takes place. 
Thus there is the constant nonzero probability density fiux j that means an possibility 
of accumulation of the wave function in close to the origin of coordinates area. In the 
second case (7 < 1/4) the collapse of the "particle" on centre misses and j = 0. 
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Further, the problem of finding of coefficient of above barrier refiection R from the 
cosmological singularity of the wave function of the Universe \E' was solved. For this 
purpose two approaches were used: 1) semiclassical approximation and 2) finding R 
as the ratio of amplitudes of wave functions of refiected and ingoing waves on infinity. 
Both approaches give identical results at 7 ^ 1/4. Thus R < 1 that means partial 
penetration of the wave function into close to zero area and its further collapse on 
centre. In case of 7 < 1/4 the second approach gives = 1 as against semiclassical 
one. It speaks about an inapplicability the latter in such situation and about necessity 
to use asymptoticses of the exact solutions. 

Let's note that on coefficient of above barrier refiection the infiuence of Hubble 
parameter has not an effect in any way because it is only renormalaze the scale factor. 

2. The model with a variable scalar field with potential V{(f) = rn^ip^ /2 was 
considered. In view of complexity of finding of an exact analytical solution the given 
problem was explored numerically. With this purpose the semiclassical solution of 
Eq. (p8D was found, therefore Eq. (pQ]) is obtained for which the system of the 
characteristics (^) was obtained. The last one was used for finding of relationship 
of the coefficient of the above barrier refiection R from value of the field $. The 
obtained relationship shows that R is rather great at Planck field. It means that 
after refiection there is an possibility of an output of the Universe on a rather long 
infiationary stage providing increase of its size with Planck up to macroscopic and a 
further output on the standard cosmological scenario. 

3. In case of nonzero factor ordering there is essentially new possibility of a 
"bounce off" a wave function of the Universe from a repelling potential wall ensured 
with a form of effective potential (|33D. By analogy to the previous case the solution of 
Eq. (P^ was found with the help of system of the characteristics (BH). The obtained 



results show that the Universe beginning the evolution with small initial value of $ 
at the stage of contraction gathers a field up to Planck. After bounce off the field 
will increase still and further Universe transfers to an infiationary stage (see Fig. 5). 

Note that in considered cases above barrier refiection and bounce of a wave func- 
tion of the Universe feature process of the quantum creation of the Universe. The 
further evolution of the model represents a stage of expansion with prompt losses of 
an anisotropy and transition into Friedmann Universe. 
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